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Abstract 

The Pauli Exclusion principle plays an essential role in the structure of the current 
universe. However, in an exactly super symmetric (susy) universe, the degeneracy of 
bosons and fermions plus the ability of fermions to convert in pairs to bosons implies 
that the effects of the Pauli principle would be largely absent. Such a universe may 
eventually occur through vacuum decay from our current positive vacuum energy uni- 
verse to the zero vacuum energy universe of exact susy. It has been shown that in such 
a susy universe ionic molecular binding does exist but homonuclear diatomic molecules 
are left unbound. In this paper we provide a first look at covalent binding in a susy 
background and compare the properties of the homonuclear bound states with those of 
the corresponding molecules in our universe. We find that covalent binding of diatomic 
molecules is very strong in an exact susy universe and the interatomic distances are in 
general much smaller than in the broken susy universe. 

PACS: 71.30 +h 



I Introduction 

String theory suggests that our current universe with a positive vacuum energy will in 
the future make transitions to other minima of the effective potential. For a review see 
ref. [1] .This prediction acquires added cogency from the essentially established fact that 
such a transition was made from an inflationary era in the very early universe to our current 
universe with a small but non-zero cosmological constant. Among the possibly numerous 
local minima of string theory |2| are some with exact susy and zero vacuum energy such as in 
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the five original superstrings. If tlie universe falls into one of these minima it might remain 
there indefinitely since there are no vacuum fluctuations in these vacua. It is therefore of 
potential physics interest to contemplate the properties of bulk susy matter. We envision a 
scenario similar to that of chain inflation [3] in which a susy ground state of the universe with 
zero vacuum energy seems plausible. However, as exemplified in this paper, in bulk matter 
the energy advantage of a system of bosons relative to a system of fermions with similar 
interactions far outweighs the vacuum energy advantage. Even apart from the possibility of 
a phase transition to exact susy, the calculation described in this paper could serve in an 
academic way to elucidate the significance of the Pauli principle in our universe. 

We assume that our current universe has a broken susy with masses of the bosonic 
quarks and leptons at the hundred GeV to one TeV scale. This may be confirmed in the 
near future at the Large Hadron Collider (LHC). Since, as in Lagrangian Higgs models, the 
number of degrees of freedom is preserved in a transition between local minima, we would 
expect that any future susy universe would have degenerate fermionic and scalar electrons, 
neutrinos, quarks, nucleons etc. In graphs and tables we assume for definiteness that the 
common mass of these susy multiplets is equal to that of the corresponding standard model 
particles in our broken susy universe where the electron has mass m = 0.511 MeV. Other 
assumptions for the common mass could be investigated since we give simple expressions for 
the Hamiltonian expectation value of the homonuclear diatomic molecules as a function of the 
common electron/selectron mass in a simple variational calculation with two parameters. 

In the following section II we present this calculation. We use a generalization of the 
Heitler-London (HL) wave function [4J first employed in the case of the hydrogen molecule, 
H2. A generalization to heavier elements is complicated in our universe due to the Pauli 
principle. However, in a susy universe where electrons freely convert in pairs to bosonic 
electrons, the generalization of the HL wave function could provide comparable accuracy for 
the binding of the heavier elements of moderate atomic number as the original calculation 
gave for In the summary and discussion presented in section HI we discuss several 
obvious improvements of the variational wave function which could be made to improve the 
accuracy of the binding energies. At present, however, it is not our intention to perform a 
state-of-the-art molecular physics calculation but to present some simple results which are 
unlikely to be qualitatively invalidated by later refinements. 

Certain integrals which occur in the calculation are given analytically in appendix A and 
a method for calculating combinatoric factors is described in appendix B. 

II Variational Wave Function for Diatomic Molecules 

For a multi-electron system of N electrons in the presence of a nucleus of charge Z the 
electronic Hamiltonian, including electron-electron repulsion, is 

H = -^T.^l-Ze^j:- + e'T.T^, ■ (2-1) 
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In a recent paper [5] a variational wave function for a susy atomic ion of atomic number 
Z with iV susy electrons was given in the form: 
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with 

uif^ = ^e-^^\^^^\ . (2.3) 

This corresponds to putting all N electrons into a Is wave function which is a poor 
approximation for > 2 in our world due to the Pauli principle. However, in a susy 
world it could be a quite adequate approximation to the ground state wave function since 
fermionic electrons in excited states would convert in pairs to bosonic electrons (selectrons) 
which would then drop into the Is ground state via photon emission. To account for the 
mutual repulsion between electrons, Zg in eq. l2.3l is treated as a variational parameter which 
minimizes the expectation value of the Hamiltonian by taking the value 

Z, = me\Z-^{N-l)) . (2.4) 

Therefore the effective nuclear charge seen by an electron is reduced from eZ to eZ^ff = eiZ— 
[N — 1)5/16) due to screening by other electrons. The radius of the atom is correspondingly 
increased by a factor of Z/Z^jf. In treating the covalent binding of two such atoms we will 
use the wave functions of eq. l2.3l centered on two nuclei at ±R. We shall see that the value 
of Zs that minimizes the energy will differ slightly from eq. l2.4l 



The Bohr radius for hydrogen is 

1 



ao = 0.529A° . (2.5) 



We work in a system of units where h = c = 1 so 

2R^ = = me^ (2.6) 

maQ 

Roo being the Rydberg constant, 13.6 eV, and being the fine structure constant 1/137. 
The variational estimate of the atomic energy is then 

E{Z,N) = -R^NZ^^ff = -R^N{Z-5{N-l)/16Y . (2.7) 

For the neutral atom (A^ = Z), Z^ff > llZ/16 oc Z. Hence with increasing Z, the atoms 
become smaller as the radius decreases as ~ 1/Z and the binding energy rapidly increases as 
Z^. As discussed in section HI, further refinements of the variational wave function are not 
expected to change these results qualitatively. From Eq. (2.7) we see that energy involved 
in subtracting or adding electrons to a neutral atom increases rapidly with Z. For example, 
the ionization energy of K{Z = 19) is 105 eV, and the electron affinity of Cl{Z = 17) is 146 



eV. In contrast, the corresponding energies in our universe are of the order a few eVs, and 
do not depend strongly Z. 

From 12.71 it is a simple matter to estimate the classical ionic binding energy by bringing 
two such atoms to a distance of twice the atomic radius, taking / electrons from one atom and 
adding them to the other taking into account the consequent Coulomb attraction. In general 
the ionic binding is much greater in the susy system than in the standard model molecules. 
We find that the binding energy increases rapidly with increasing AZ = |Zi — Z2I. This 
is because it is energetically favorable for a low-Z atom to give up an electron or two to a 
high-Z atom since ionization energy and electron affinities increase rapidly with Z . However, 
the ionic binding energy is found to be negative (anti-bonding) for homonuclear molecules 
(AZ = 0) in both the susy and standard model systems. 

Covalent bonding is an intrinsically quantum mechanical effect where one or more elec- 
trons are simultaneously shared by two nuclei. In a susy system the lowest energy will be 
found by the totally symmetric wave function with at most two fermionic electrons. 

Our approach is to follow the early quantum mechanical treatments of molecular bonding 
while discounting the effects of the Pauli principle. We consider a diatomic system of two 
atoms each having nuclear charge Z and N = Z electrons. 

One might try the many electron generalization of the molecular orbital method [6] 
corresponding to the wave function 



Here 2R = D is the interatomic distance to be determined variationally. Although this gives 
a rough approximation to the H2 molecule it would not be a serious candidate for a many- 
electron molecule in our universe due to the Pauli principle. In the susy case it becomes of 
interest to consider this wave function for larger A^. The wave function 12.81 will be an exact 
solution in the absence of e-e interaction. However, for small AZ it does not work well since 
the wave function, when multiplied out, contains a preponderance of ionic configurations in 
which one atom has more electrons than the other. These are not favorable energetically 
because of electron-electron repulsion within the atom, and as noted above, lead to unstable 
ionic bonding for homonuclear molecules. Further, since the weights of these configurations 
do not depend on R, the total energy in the wave function of eq. l2.8l does not approach 
twice the energy of the isolated atoms given by eq. l2.7[ In the actual state the electrons' 
motion will be correlated to keep them apart. Indeed, for moderate values of Z, the energy 
expectation value for the wave function eq l2.8l is found to be greater than that from the 
generalization given below of the Heitler-London wave function [1] which therefore provides 
a better approximation to the true wave function. 

We consider, therefore, avariational wave function of the form: 



2N 




(2.8) 



i=l 



vl> = No{N) Yl u{fi - i?) Yi "(^i + R)\ + permutations {fi ^ fj) . (2.9) 



Here Nq{N) is a normalization constant. Note that there are equal number of electrons 
centered on each atom. The permutations make the wave function totally symmetric as 
appropriate to bosonic electrons. If there are (up to two) fermionic electrons in the system, 
they would also be expected to have a ground state in the totally symmetric state since 
other possibilities would be partially anti-bonding. We do not discuss excited states in the 
current paper nor do we consider possible mixing effects between states of differing number 
of fermionic electrons. In the presence of the other atom, Zg would be R dependent and 
would be expected to approach that of the isolated atoms given by eq. l2.4l as R oo. We 
treat Zg and R as variational parameters to minimize the energy. 

In the general case, the symmetrization of the wave function leads to a complex combina- 
torial problem involving a large number of terms each of which is a product of 2N m's. These 
can be grouped into configurations with N + m electrons in the wave function, u = u{\f—R\), 
centered on the nucleus at R and N — m electrons centered at —R. The solution of this 
problem is given in Appendix B. Here we specialize to m = which corresponds to the wave 
function of eq. l2.9[ In this case there are (2A^)!/A^!^ permutations, including the stated refer- 
ence configuration where the first N electrons are centered at R and the second N electrons 
are centered at —R. Among the permutations, the number of terms with k interchanges is 
the square of the binomial coefficient 

n(A;)=(^^y . (2.10) 

Note that 

N 

J2n{k) = {2N)\/N\^. (2.11) 

k=0 

As a function of Zg and R we minimize the expectation value of the Hamiltonian written as 
a sum of kinetic, electron -nucleus Coulomb terms, electron-electron correlation terms, and 
nucleus-nucleus repulsion term. 

1 2N 

Hkin= (2.12) 

2N 1 

r,- — rA 

i<j I * J I 

H^uc=Z^ey\2R\ (2.15) 

The corresponding expectation values are written in terms of one-body and two-body in- 
tegrals, Ai, summarized in appendix A together with combinatoric factors discussed in ap- 
pendix B. Each expectation value takes the form of a sum over k, the number of interchanges. 
The normalization of the wave function, for example, requires that 

(NoiN))' = (E^(^)^o') (2.16) 



The results are 



< mk^n\'^ >= iNoiN)f — E ^(k) ({N - k)Af + kA,Af~') (2.17) 

\ / fc=0 



(2.18) 



< vl/|if,,|v[/ >=(iVo(iV))' [^!!\ Z,e^ J2 <k) [A;,Af~^k{2k - 1) + Ak{N - k)A^,Af'^ 

+ Af{N - k)ih{N -k-l) + A,,{N - k))) . (2.19) 
In the large R limit 

< <i/\H\^ >^ N (Zl/m - 2e^Zs{Z - b{N - 1)/16)) (2.20) 
This is a minimum at 

Zs = me^[Z - 5(Ar - 1)/16) = {Z - 5{N - l)/16)/ao (2.21) 

with 

< H -2R^N{Z -5{N -l)/lQf . (2.22) 

Thus at large separation the energy of the diatomic system approaches twice the energy of 
each isolated atom given in eq. l2.7[ By examining the asymptotic forms for the integrals Ai 
one can see that for large R the asymptotic energy is generally approached exponentially. 

At intermediate values of the separation we write 

Zs = (Z -5{N -l)/16){l + 6)/ao . (2.23) 

We vary 6 and the interatomic separation 2R to minimize the energy. In figure 1 we 
show for H2 the expectation value < > as a function of R at the optimum value of 6 
for each R. In table 1 we give the results of the overall minimization for the eight lowest 
homonuclear diatomic molecules. We show for each molecule the minimum energy and 
equilibrium interatomic separation as well as the overall optimum value of 6. We neglect the 
vibrational energy of the nuclei in the ground state so we do not distinguish between the 
minimum value of the Hamiltonian and the negative of the dissociation energy. 
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Figure 1: For the case of the O2 molecule with bosonic constituents the graph plots the 
expectation value of the Hamiltonian in eV as a function of the interatomic distance D in 
A°. For each value of D, the 6 parameter is chosen to minimize the energy. 
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H2 


0.748 


-3.782 


0.166 


0.74 


4.71 


He2 


0.454 


-27.15 


0.208 


2.97 


0.0009 


Li2 


0.326 


-86.04 


0.222 


2.67 


5.4 


Be2 


0.254 


-196.3 


0.228 






B2 


0.208 


-374.0 


0.231 


1.59 


3.0 


C2 


0.176 


-635.0 


0.234 


1.24 


6.21 


N2 


0.153 


-995.6 


0.235 


1.10 


9.76 


O2 


0.135 


-1471.6 


0.237 


1.21 


5.12 



Table 1: Binding properties of the eight lightest homonuclear diatomic molecules with 
totally symmetric electronic wave function as would be appropriate in the susy case. In the 
second column we give the interatomic distance {D = 2R) in angstroms. The corresponding 
minimum energy expectation value is given in eV in the third column. The fourth column 
gives the 6 parameter at the minimum. The fifth and sixth columns give the experimental 
[7] interatomic distance in angstroms and the experimental binding energy in eV for the 
diatomic molecules of the current universe (fermionic constituents). Neglecting vibrational 
energy the experimental binding energy is the negative of the potential at the minimum. 



Ill Discussion 



We have studied the molecular properties of the first eight homonuclear diatomic molecules 
in a totally symmetric HL electronic ground state as would be appropriate in a susy back- 
ground where fermionic and bosonic electrons are degenerate in mass. This wave function 
was used for H2 by Wang [S] in the early days of quantum mechanics. 

Our results for H2 are in good agreement with those of that paper obtained in the days 
before computers. For instance our results from table 1 for the H2 case can be compared 
with his interatomic distance of 0.76 A° , minimum energy of —3.76 eV, and 5 parameter 
0.166. The H2 calculation, but not that for heavier molecules, is the same for bosonic as for 
fermionic electrons due to the absence of Pauli effects. 

Our model calculation and that of Wang is 20% below the experimental binding energy 
of 4.71 eV. For the interatomic distance these calculations are good to about 1%. We can 
expect that similar accuracies would hold in the susy case for the diatomic molecules above 
H2- The interatomic distances predicted in table [T] are somewhat less than twice the radii 
of the isolated atoms ?>/{2Zs). 

In the case of the higher molecules, our calculation shows progressively smaller in- 
teratomic distances compared to the experimental values for the molecules with fermionic 
electrons which are all of order of 1^1°. The third column of table 1 shows the energy of di- 
atomic molecule relative to the total energy of two isolated atoms. We see that molecules are 
much more strongly bound than the corresponding molecules with fermionic constituents. 
The binding energy increases rapidly with Z. In contrast, in the standard model molecular 
binding energy is of the order of few eVs in all cases (column 6). These facts suggest that 
molecules in a susy background would have much lower reaction cross sections than the 
diatomic molecules in our broken susy background. 

From the second derivative of the energy at the minimum one can estimate the first 
vibrational energy level. 



where M is the nuclear mass. From the curve of figure [T] we can estimate that the vibrational 
ground state lies, in the case of O2, 4.4 eV or 0.3% above the minimum of the potential. 
This compares with 0.196 eV in the case of the usual O2 molecule. [7]. 

In summary, the tendency of bosonic electrons to accumulate in the low-energy regions 
leads to strongly bound atoms and molecules in a susy world. In a susy atom all electrons 
occupy the Is orbital in the ground state, leading to a spherical shape whose radius decreases 
as l/Zeffi and the total energy scales as - ^e//- The important point is that nuclear charge 
is only partially shielded, so that Ze// oc Z, which leads to ionization energies and electron 
affinities that increase rapidly with increasing Z. Similarly, atoms strongly bind to form 
molecules with binding energies that increase rapidly with Z. For example, a binding energy 
of 1471 eV for O2 corresponds to a dissociation temperature of ~ lO^i^'. These properties 
have important consequences for bulk matter in the susy universe. It suggests that com- 
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plex molecules with larger number of atoms would form easily and that solids might have 
high melting points. Of course, macroscopic systems would exhibit Bose condensation and 
superfludity with high critical temperature. 

The situation is very different in our universe, where only valence electrons take part 
in chemical reactions and molecular bonding, as they occupy the outermost orbitals due to 
the Pauli principle. The other electrons are effectively frozen out as they occupy the inner 
orbitals; their primary function is to shield the nuclear charge, which they do quite efficently. 
Consequenly Zg/j ~ 1, and ionization energies and electron affinities, and molecular binding 
energies are all of order a few eVs for all atoms. 

As in every variational calculation, the accuracy can be improved and the energy estimate 
lowered by adding additional variational parameters. Asymptotically, the exact ground state 
energy will be approached from above. In the case of H2 an improved energy was obtained [H] 
in the early days by adding a single parameter corresponding to a term in the wave function 
where both electrons were centered on the same nucleus. Reference [9j found a value of 4.0 
eV for the dissociation energy compared to the experimental value of 4.7 eV. 

For higher molecules similar refinements of the wave function would lead to an improve- 
ment in energy, but the main results are not expected to change qualitatively. In the case 
of a single atom, the key issue is the screening of the nuclear charge. Obviously the result 
can be improved by using the self-consistent Hartree-Fock approximation, in which the one- 
particle wave function u{r) in 12.21 itself is determined variationally. Because of the spherical 
symmetry, the screening electric field seen by an electron at R is due to the fraction of the 
other N — 1 electrons that are within the sphere of radius R. Hence the average screening 
charge is given by 



Since the function u{r) is the same for all electrons, a finite fraction of each electron is outside 
the sphere. Only when R ^ 00, does the integral equal unity, and Qscr = —{N — l)e. On the 
other hand, suppose we take R to be the most probable distance, obtained by minimizing 
r'^v?{r). Then roughly half the particles will be outside the sphere of radius R, and hence 
roughly half of Z nuclear charge will be unscreened. This means that, Ze// oc Z. Note 
that this argument does not depend on the precise form of u{r). Hence, our results are 
qualitatively correct. Similar quantitative improvement can be achieved by introducing two- 
body (so called) Jastrow factors of the form Hij /(^i ~ "O) wave function, which take 
into account correlation effects by keeping particles away from each other. The problem 
then becomes intractable for large Z. However, we do not expect the qualitative physics to 
change. 

For molecules results can also be improved by including admixtures of ionic configurations 
- that is, configurations with Z+m particles centered on one atom, and Z—m atoms centered 
on the other, with the relative amplitude determined variationally. These effects will lower 
the energy, and thus make the molecules even more strongly bound. However, since the 
transition to exact susy is unlikely to occur in the near cosmological future, it does not seem 
urgent at this time to seek an accuracy beyond that obtained by our present approximation. 
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A Overlap Integrals 

Since the wave function is totally symmetric in the 2N electron positions, overlap integrals 
can be written as the number of permutations (2A^)!/A^!^ times the overlap of a reference 
configuration with the full wave function. For the reference configuration we take the wave 
function where the first (scalar) electrons are in Is states concentrated around the nucleus 
at position R while the remaining electrons are concentrated at —R. 



m^^f = No{N)(\{u{n-R)\ [ n u(r; + i?)] =!--...- + + ...+ > (A.l) 

\i=l / \j=N+l ] 

Then the matrix element of any operator takes the form 

<%0\^>=NQ{Nf <^^^^\0\^^^^ Sperms > . (A.2) 

If (9 = 1 we have the normalization condition of eq. l2. 161 with 

=< -1+ >= S / d^r^-Zslr-Rl-Z^r+Rl ^ ^-2R' ^ ^R' + AR'^S) (A.3) 
71 J 

where 

R' = ZsR (A.4) 

The Hamiltonian consists of a sum over single particle terms and two particle terms. Its 
expectation value is given by a combinatoric sum over single and double particle expectation 
values. For instance 

■e-2H(l-2/|r1) = 1 (A.5) 
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and 

=e-2^'(l + 2i?'-4i?'V3) (A.6) 

The electron-nucleus interaction involves the integrals 

h=< -| ^ ^ ^ |- > =< +1 ^ ^ ^ 1+ >= 1 (A.7) 
\r — R\ \f+R\ 

Ma=< - L ^ ^, 1- >=< + L ^ ^, 1+ > (A.8) 

\r + -R| \r — R\ 

A2b=2<-| — ^|+>=2<+| — ^|-> (A.9) 
f±i? \r±R\ 



-I r\ 



It is easy to show that 

^2a=-l + (l + ^)(l-e-^^') (A.IO) 
A2fe=2e-2^'(l + 2i?') (A.ll) 

The electron-electron repulsion term involves the two body integrals 

/3=< - - 1 1 ^ ^ ^ 1 1 - - >=< + + ^ ^i l + +> (A.12) 
|ri - r2| |ri - r2| 

A^a=< - + 1 1 ^ ^ ^ 1 1 - + >=< + ^ ^i l + -> (A.13) 

|ri-r2| ^1-^^21 

Azh=< - + 1 1 ^ ^ ^ 1 1 + + >=< - + | |^ ^ ^| |--> (A.14) 
|ri-r2| ri-?"2| 

A^^=< - - 1 1 ^ ^ ^ 1 1 + + >=< + + ^ ^| |--> (A.15) 
|ri-r2| In -^2! 

All of the integrals of this appendix except for occurred in the analysis of the H2 molecule 
in the 1920's and were reviewed in the book of Pauling and Wilson [TD]. We have added the 
evaluation of A'^^ since it occurs in the binding of the higher diatomic molecules in the susy 
case. The analytic expressions are 

/3=5/8 (A.16) 
^3a=^(l - e-'^') - e-^^'(ll/8 + 3i?72 + 2B!^ /?,) (A.17) 

A,,=e-'^' (^^^(5 + 4i?') + 2i?'j (A.18) 

and 

A3,=-e-^^' (25/8 - 23i?V2 - 12i?'' - 8i?'V3 + + 2B! + AB!^ I?,f{-i + ln(2i?')) 

+ _|.(e8^'^,(_8i?')(l - 2B! + 4i?'V3)2 

-2e''^' Ei{-AB!){1 - 2B! + 4i?'V3)(l + 2B! + AW^ /?,))) (A.19) 

The exponential integral function Ei(x) goes asymptotically as which makes it convenient 
to factor out the overall e~^^ in eq. IA.19l The complete expression for 743c was given in 1927 
by Sugiura pT] . 

B Combinatorics 



We consider a diatomic molecule with the two atoms centered at R and — i?, each with 
nuclear charge Z, with a total of 2A^ = 2Z electrons. However here, for greater generality, 
we take Z and to be independent. The many-body wave function ^ is constructed from 



single-paricle functions ^(r ± R) centered at —R and R, respectively. The function u{r) is 
similar, but not identical, to the atomic wave function, as it is calculated variationally. 

The total wave function \E' is a linear combination of functions ipm having N + m particles 
centered at —R and N ~ m particles centered at i?, with — < m < N: 



N-m _^ 2N 

1pm = Y\. '^i^i ~ n uifj + R) + permutations . (B.l) 

i=l j=N-m+l 

In calculating the expectation values we need to solve a combinatorial problem. This can 
be done in a compact way by introducing a pseudospin variable /i for each particle such that 
;U = 1 or — 1 if the particle is centered at —R and -R, respectively. Then Y^^iHi = 2m. 

Let Uf^{r) = u{f+ fxR). 

Then the ground-state wave function is given by 

/2jY\ -1/2 TV 2N 



\ / m=-N {/i} ^ i=l 



Where the a^'s are variational parameters. By symmetry am = O-m- 

Consider first the evaluation of the normalization constant < \E'|\E' > which has 2^^ terms, 
each of which is an integral over the 2N coordinates, and thus consists of 2N factors. We 
assume that each function M^(r) is normalized. But, in general, and has a non-zero 
overlap Aq, so that 

j £ru^{r)u^>{r) = 5^^^> + AqS^,.^/ . (B.3) 
Then the normalization constant has the general form 

2N 

<^|^>=E^p^o (B.4) 

where Zp is a combinatorial factor. To determine Zp consider orthonormal functions 

t;n(g)=e2™« (B.5) 
where < g < 1, and n is an integer. Since 

< VmK >= C dqe^^'^^-^^^ = 5n,m , (B.6) 

Jo 

we can represent the Kronecker 5 in the wave function IB. 21 in a similar manner, leading to 



The representation allows us to integrate over and sum over /Xj for each i separately, which 
gives 



V / m,m' ^ ^ 



Expanding the factor {...}^^, and integrating over q,q' and after some rearrangements we 
finally obtain, 



2N-\m+m'\ 

<*|^>=^a^a^, ^ F{N,p,m,m')Al . (B.9) 

m,m' p=|m— m'l 



where. 



with 

G = {A^ + (m + m' - p)/2}!{iV - (m + m' + p)/2}!{(p + m - m')/2}!{(p + m' - m)/2}\ . 

Thus F is zero unless p + m + m' is even. 

Ground-State Energy: To evaluate the energy < >, let us write H — Hi + Vee, 

where Vee describes the electron-electron interaction, and 

Hi = ^21=1 Hi is the sum of one-particle Hamiltonians Hi. We first consider Hi whose 
average is given by < \E'|ifi|\E' >= 2N < \E'|iJj|^ > since the particles are identical. The 
evaluation of < > is similar to that of < > since the integration over the 

coordinates factorizes. Let 

j SriU^{ri)HiU^i{ri) ^ Ai5^t,> + A25^^-fj,> . (B.ll) 

Proceeding as before, we obtain, 

< m\Hi\^ >= amam'I(m,m'), (B.12) 



where 



2N-\m+m'\ 

I{m,m')= F{N,p,m,m'){{2N-p)AiAPo+pA2A^o~^} . (B.13) 

p=\m—m'\ 



Then the expectation value of the electron-electron term is given by 

since all pairs yield the same expectation value. The procedure is similar, except that 
averaging over the pair interaction involves integrating over two coordinates. Let 



B^= f dh^Sr2ul{n)-^^-^ul{f2) . (B.15) 

I ^3 I 

B2= f d^n£r2ul{n) ^ | H-^(r2) . (B.16) 

I ^3 I 

2 

B^= j d^nd\2ul{fi) \ ^ u^{r2)u.^{r2) . (B.17) 

I ^3 I 

2 

S4=/dVidV2ii^(fi)li_^(fi)— ^M^(f2)«_^(f2) . (B.18) 



Then we have 



< *|i;,e|* >=l^l + l^2 + ^3 + ^4, (B.19) 

where 

1N-\m-\-m'\ 

Vi^ F{N,p,m,m'){{2N -p){2N -p-2)^{m^m'f}BiAl/A . (B.20) 

p=\m—m'\ 
2N-\m+rn'\ 

J2 F{N,p,m,m'){{2N -pf - {m + m'f}B2Al/A . (B.21) 

p=\m—m'\ 
2N-\m+m'\ 

1^3= E i^(A^,P,m,m')(2iV-p)p53Ag-' . (B.22) 

p=\m—m'\ 
2N-\m+m'\ 

V^= F{N,p,m,m')-^{{^p{p-2) + {m-m'f)Al-'^ 

p=\m—m'\ 

+ ((2N -pf -{m + m'f)Al} . (B.23) 

Heitler-London limit: In this case we keep only m = m' = 0, so that each atom 
contains exactly N electrons. As i? — > oo, the molecule dissociates into two isolated atoms 
in their respective ground states. We can set oq = 1. 

Then we have 

<^|^>=^ Af. (B.24) 

fe=0 V ^ / 



< E , {2(^ - k)AfAi + 2kAl''-^A2}. (B.25) 

fc=0 V ^" / 

< *|ye,e|* >=E ^ J (Ag'{(A^ -k){N-k- + (TV - A;)2(S2 + B,)} 

+ 4k{N - k)Af-^B^ + k{k - l)BiAf-^) . (B.26) 

If we include m,m' = ±1, there is one extra variational parameter Oi = a_i. These 
contribute the following terms to < \E'|\E' >: 

This analysis is given for completeness and as a guide to future improvements only. In the cur- 
rent paper, numerical results are presented for the m — m' — term only. 
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